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Yang-Mills color fields evolve chaotically in an anisotropically expanding universe. The chaotic behaviour differs from that 
found in anisotropic Mixmaster universes. The universe isotropizes at late times, approaching the mean expansion rate of 
a radiation-dominated universe. However, small chaotic oscillations of the shear and color stresses continue indefinitely. An 
invariant, coordinate-independent characterisation of the chaos is provided by means of fractal basin boundaries. 


Yang-Mills fields are central to quantum theories of el¬ 
ementary particles. They are of interest to dynamicists 
since they evolve chaotically in flat spacetime 
But how do Yang-Mills fields behave in the early uni¬ 
verse? Does the chaos persist, or is it eradicated by 
the general relativistic effects of cosmological expansion? 
The earliest studies of the Einstein-Yang-Mills (EYM) 
system assumed isotropic cosmological expansion. This 
imposes special symmetries on the dynamics and the dy¬ 
namical system is integrable: chaos cannot exist Q. Re¬ 
cently, a Yang-Mills theory was formulated in an axisym- 
metric, spatially homogeneous universe |^. The large 
number of degrees of freedom made an analysis of the 
full dynamics difficult, and the coordinate dependence of 
the standard chaotic indicators meant that the relativis¬ 
tic chaos could not be invariantly characterised. 

In this letter, we analyse general relativistic Yang-Mills 
chaos using invariant topological methods. We make the 
problem more tractable by an economical definition of 
variables, which reduces the dimension of the phase space 
substantially (from 8 — D to 5 — D). In the physical pic¬ 
ture that emerges, the asymptotic evolution of the spatial 
volume of the universe imitates a radiation-dominated 
universe, while the shear diminishes chaotically. 

Previous studies of chaos in general relativity have fo¬ 
cused on the non-axisymmetric Bianchi type VIII and IX 
(Mixmaster) universes, where the presence of anisotropic 
3-curvature creates an infinite sequence of chaotic oscil¬ 
lations on approach to an initial Weyl curvature singu¬ 
larity at t = 0 This behaviour is intrinsically 

general relativistic. By contrast, the chaotic EYM cos¬ 
mology that we study is different: chaos exists even when 
the metric is axisymmetric and the curvature is isotropic. 

We evolve the Yang-Mills fields in the simplest 
anisotropic metrics of Bianchi type I. The color degrees 
of freedom of the Yang-Mills gauge fields oscillate chaoti¬ 
cally, while the expansion attenuates their overall energy. 
The EYM action is 
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where F^i, is the gauge-invariant field strength. The 
spacetime is described by the axisymmetric Bianchi I 
metric, with scale factors b{t) and c{t): 

ds^ = —dt^ h^{t) {dx^ -\- dy^) -b (?{t)dz^. (1.2) 


After fixing the internal gauge of the Yang-Mills held, 
the matter can be parametrized by two variables (a, 7) 
Q, which can be thought of as color degrees of freedom 
for the massless gauge fields. The cosmological evolution 
is an orbit in the (6, c, a, 7) phase space. 

We define the mean expansion scale factor by a = 
(6^c)^/^, the shear anisotropy by % = {b/c)^/^, with 
volume expansion Ha = d/a, and shear = x/x- 
The scaled Yang-Mills field strengths will be defined 
by (4',r) = (a/(ax), 7x^/a), with conjugate momenta 
(n3i,nr) = (Q!/(ax),7x^/a). The Einstein equations re¬ 
duce to 


Ha+2H^ + H^ = 0 (1.3) 

H^-Hl + Hl + 3HaH^ = ^A^^ + ^Ul . (1.4) 

The Hamiltonian constraint equation is 

+ + = 3{Hl-Hi) . (1.5) 

The conservation of gravitational-plus-matter energy ap¬ 
pears in the (4>, T) subsystem as a loss of energy, given by 
E{t) = 3{Ha — Hi). The matter conservation equations 
are 

n,i, -b 2{Ha + H^)U^ -b T- (r^ -b 4-2) = 0 

ilr + 2{Ha + H^)Ur + 2A>‘^T = 0 . (1.6) 

The matter sector is thus a driven, dissipative system. 
The Yang-Mills coordinates decay adiabatically due to 
the expansion, but they are also driven by oscillations in 

Hx- 

There are 4 coordinates, {Ha, Hy,, 4*, T), and 8 degrees 
of freedom. This system permits at m ost 4 positive Lya¬ 
punov exponents. However, eqns ( 0 ) and (0 are con¬ 
straints. Consequently, only 6 constants of integration 
need be specified. The Hamiltonian constraint reduces 
this by one, leaving 5 degrees of freedom. If each of 
the three constraint equations are unique, then there is 
only one positive Lyapunov exponent, as in flat space- 
time {Ha = H^ = 0). This is to be expected since the 
vacuum Bianchi I universe is integrable. Any chaos in 
the metric variables is therefore a consequence of chaos 
in the Yang-Mills field. This was also found to be the 
case in Ref. |§]. 


1 





The color amplitudes of the Yang-Mills system scatter 
around the potential of Fig. |^. The chaos is more com¬ 
plicated than that found in Mixmaster universes because 
of the rapidly varying curvature of the hyperbolic poten¬ 
tial walls. In flat spacetime the energy remains constant, 
and the contours define different surfaces of constant en¬ 
ergy. In the Bianchi I universe o f ge neral relativity, E{t) 
drops asymptotically (see eqn (1^)). The four hyper¬ 
bolic potential walls shrink as the energy available to the 
Yang-Mills field decays. We see from Fig. ||, the trajec¬ 
tories occupy smaller and smaller volumes of the matter 
phase-space. 

Chaos is often quantified by computing Lyapunov ex¬ 
ponents. However, the values of Lyapunov exp^onents are 
coordinate dependent in general relativity 0, because 
of the coordinate covariance of Einstein’s equations (an¬ 
other manifestation of the so called ‘problem of time’ 
in cosmology). The authors of Ref. 1^ cite this fact, 
together with the large number of degrees of freedom 
(which lead to Arnold diffusion) and the non-compact 
phase space, as major obstacles to a generalization of 
the study of chaotic properties of Yang-Mills theory to 
curved spacetimes. These barriers can be overcome by 
the methods of chaotic scattering which al¬ 

low us to identify fractal sets which fully characterize the 
chaos. Fractals cannot be hidden by a reshuffling of co¬ 
ordinates: their existence and dimension are coordinate- 
independent topological features ||,^,|^,|lT|. The frac¬ 
tal set that we seek is the set of all periodic orbits: by 
analogy with the accessible states in thermodynamics, 
this periodic set completely describes the chaotic dy¬ 
namics. The fractal set is also called a ‘strange repel- 
lor’, or ‘strange saddle’. If the global expansion of the 
universe could be projected out of the dynamical evo¬ 
lution, the repellor would consist of all periodic orbits 
trapped in the potential. When the global expansion 
is included, the periodic orbits are actually self-similar, 
reminiscent of the dynamics of the Bianchi IX cosmology 


The strange repellor can be coded by a symbolic dy¬ 
namics, or isolated directly by numerical techniques such 
as the PIM method. It is easier to employ the method 
of fractal basin boundaries. A slice of initial conditions 
is taken through phase space. All possible asymptotic 
states are determined and assigned a color (black or white 
here). If large blocks of initial data space suffer the 
same fate, then the outcome basin will look smooth and 
monochromatic; by contrast, highly mixed, fractalized 
basins indicate a sensitivity to initial conditions, as well 
as mixing and folding of trajectories. Hence, fractalized 
basins signal chaos in a covariant way. All observers agree 
upon the occurrence of the events used to construct the 
fractal, and all will agree on the dimension of the basin 
boundary. 

A typical trajectory will travel down one of the F chan¬ 
nels before rebounding back into the scattering region of 
the potential. This is again reminiscent of the Mixmas¬ 
ter system, where the repellor was also inefficient 


In Mixmaster, the repelling set can be isolated by artifi¬ 
cially opening the exit pockets, so that orbits thrown far 
from the scattering region were allowed to escape. Sim¬ 
ilarly, here, we cut holes in the pockets and assign the 
color white to the initial condition if the orbit falls down 
the upper pocket (F > 0), and black if it falls down the 
lower pocket (F < 0). In flat spacetime, this is straight¬ 
forward. In the Bianchi I spacetime, the potential walls 
move inward as the energy is not constant. However, 
the same procedure can be followed, so that the angular 
size of the hole in the pocket is the same for all energy 
contours. Basin boundaries, sliced through ('I',F) coor¬ 
dinates, are shown in Fig. The pockets are of size 
6 = 2arctan0.5. The fractal basins look very similar to 
those we find in a flat spacetime. The box-counting di¬ 
mension was estimated to be maximal, namely. Do = 2 
(the same dimension as we found for the basins in flat 
spacetime). If the pockets are made infinitesimally thin, 
the fractal boundary fills all of phase space. 

The interweaving of outcomes, demonstrated by the 
basin boundary structure, corresponds to a final-state 
sensitivity that is directly related to the fractal dimen¬ 
sion. An e uncertainty in initial values leads to a final- 
state uncertainty of where N is the dimension of 

the phase space. In the N = 2 slice of Fig. ||, we found 
Do = 2. Axisymmetric, Bianchi type I, EYM cosmolo¬ 
gies are therefore very sensitive to initial conditions, and 
the shear evolution is highly chaotic. 

The chaotic scattering has a subtle effect on the large- 
scale structure of the spacetime. In vacuum, the scale 
factor grows as a cx and the shear is of compara¬ 
ble importance, x ^ The color oscillations scat¬ 

ter the metric variables chaotically. As the universe ex¬ 
pands, the volume of the ('I',F) phase space redshifts as 
p ^ 1/a^. The scale factor eventually evolves towards the 
behaviour of a radiation-dominated universe, a ^ 
and is unaffected by the chaotic color oscillations. The 
relaxation timescale for this behaviour is related to the 
Lyapunov exponent, and hence to the fractal dimension 
[p^ . A rough estimate indicates that trajectories at¬ 
tract onto Ha = l/(2t) after a few e-folds of the scale 
factor. The shear is more vulnerable to stochastic be¬ 
haviour and the color oscillations sustain oscillations of 
about zero. This general behaviour can be seen by 
examining the scattering angle in minisuperspace. A tra¬ 
jectory leaves the (lna,lnx) minisuperspace plane at an 
angle 0 = a,rcta,ii{H^/Ha). Fig. |4| compares the scatter¬ 
ing angle in the case of regular motion to the angle in 
the chaotic case. The scattering angle is not constant in 
time as can be seen in Fig. While Ha escapes the ef¬ 
fects of the chaotic scattering, H^ remains sensitive and 
continues to oscillate. 

The dependence of the final scattered angle on the im¬ 
pact parameter is a direct probe of the strange repellor. 
In a simple chaotic system, the scattering angle can be¬ 
come fractal. No matter how small a difference in the 
incident velocity, the difference in scattered angle is size¬ 
able. The scattering angle and the basin boundaries are 
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different ways of viewing the same phenomenon. Regard¬ 
less, the scattering angle provides an important perspec¬ 
tive. When the motion is regular, the angle approaches 
±7r/4 and so ^ ^Ha- The evolution of the shear is 
as important as the global expansion of the spacetime 
volume, and initially the spacetime is highly anisotropic. 
At late times, under the influence of chaos, the scatter¬ 
ing angle clusters about zero, and \H^\ Ha, with the 
shear frozen at a nearly constant value. Since the field 
equations were shown to be invariant under a rescaling 
of the shear, we can choose that scale so that Inx —s- 0 at 
freeze out. Physically, this is equivalent to the statement 
that a b c. The anisotropy decays and the final state 
looks isotropic. 


The chaotic behaviour influences the solution only 
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FIG. 1. Isocontours of the color potential (eqn. (1.6)). As the energy, E{t), decreases, the walls of the potential move inward 
and the color oscillations are confined accordingly. A typical trajectory is shown. 
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FIG. 2. A projection of the Bianchi I phase space onto the ('F,n.ii) plane (left) and onto the (r,nr) plane (right). The 
initial values for this trajectory are = 0.1, F = 0, Ha = 0.11, = 0.08, ^' = 0 and F is fixed by the Hamiltonian constraint. 
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FIG. 3. The strange repellor is revealed in these fractal basin boundaries. Holes are cut in the F channels to allow typical 
trajectories to escape. The angular size of the pockets in this simulation is tan(0/2) = 0.5. 
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FIG. 4. The scattering angle defined as 0 — aTcsm{H^/Ha) as a function of the initial Ha- The left plot is not chaotic with 
Ik = Ik = 0. In the right panel the initial conditions for the remaining degrees of freedom are H^o = 0.0999,4^ = 0.1 = —F = —4/. 
Again, the Hamiltonian constraint fixes F. 
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